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Low-Energy Charge-Density Excitations in MgB 2 : Striking Interplay between 
Single-Particle and Collective Behavior for Large Momenta 
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A sharp feature in the charge-density excitation spectra of single-crystal MgB2, displaying a 
remarkable cosine-like, periodic energy dispersion with momentum transfer (q) along the c*-axis, 
has been observed for the first time by high-resolution non-resonant inelastic x-ray scattering (NIXS). 
Time-dependent density-functional theory calculations show that the physics underlying the NIXS 
data is strong coupling between single-particle and collective degrees of freedom, mediated by large 
crystal local-field effects. As a result, the small-g collective mode residing in the single-particle 
excitation gap of the B n bands reappears periodically in higher Brillouin zones. The NIXS data 
thus embody a novel signature of the layered electronic structure of MgB2 . 

PACS numbers: 71.45.Gm, 73.43.Lp, 74.70.Ad, 78.90.+t 



The discovery of superconductivity in MgB2 at a high 
critical temperature T c ~ 39 K 1] has stimulated nu- 
merous theoretical and experimental studies worldwide. 
Unlike high-T c cuprates, the relatively simple crystal 
structure of MgB 2 allows detailed first-principles cal- 
culations to be performed and compared with experi- 
ments. As a result, MgB>2 is now widely accepted as 
a phonon-mediated conventional superconductor based 
on the anisotropic Eliashberg formalism with multiple 
gaps facilitated by the different electron-phonon coirpling 
strengths of the boron a and 7r bands 0, 0, 0,0, IS 0- 
The strong electron-phonon coupling between the 2D a 
bands and the in-plane vibration of the B layers (the E^ g 
phonons) dominates the superconducting properties and 
is largely responsible for the unusually high T c . 

Within this phonon-mediated picture of supercon- 
ductivity, dynamically screened electron-electron and 
electron-ion interactions play an important role 0|. 
The study of charge-density response function through 
inelastic x-ray 0] or electron |ld( scattering experi- 
ments and/or theoretical calculations provides a power- 
ful means for investigating the details of the dielectric 
screening and the associated crystal potential, local-field, 
and exchange-correlation effects. Recent first-principles 
calculations f° r example, have shown the impor- 

tance of the crystal local-field effects (CLFE) in shaping 
the superconducting properties of MgB 2 . Earlier theoret- 
ical studies of the collective charge excitations in MgB2 
based on first-principles calculations by Ku et al. [l3j 
and Zhukov et al. 14], on the other hand, predicted a 
sharp new collective mode between 2-5 eV for momentum 



transfer g||c*-axis as a result of coherent charge fluctu- 
ations between the Mg and B layers due to the unique 
electronic structure of MgB 2 . The feature as predicted 
dispersed weakly with q and decayed via Landau damp- 
ing into the single-particle continuum at q ~ 6 nm" 1 in 
the first Brillouin zone (BZ). The CLFE were reported 
to be negligible, for the small g-range investigated. 

In this Letter we demonstrate, using state-of-the-art 
high-resolution non-resonant inelastic x-ray scattering 
(NIXS) experiments and ah initio time-dependent den- 
sity functional theory (TDDFT) calculations, that the 
charge response of MgB2 is truly remarkable. Indeed, 
the long-lived, low-energy collective excitation in MgBo 
exists not only for small q as previously predicted (l3lll4| . 
but it actually extends to higher BZ's along the c*-axis. 
The mechanism behind such physics is the strong cou- 
pling between the single-particle and collective excitation 
channels which is mediated by large crystal local-field ef- 
fects (CLFE) due to charge inhomogeneity normal to the 
Mg and B layers. As a result, the conventional Landau- 
damping mechanism does not restrict the collective mode 
to a small fraction of the momentum space: the mode ac- 
tually reappears periodically as q steps through successive 
BZ's. The enhanced scattering cross section at high q and 
the negligible multiple-scattering effects of the x-ray mea- 
surements play a decisive role for the clean observation 
of this mode over a large momentum range. Our work 
represents also the first observation of collective charge 
excitations of this kind in condensed matter physics [15J . 

The MgB 2 single crystals were prepared by high- 
pressure sintering of MgB2 powder |l£|. Superconduct- 



2 



ing transition with the onset temperature of 39 K was 
confirmed by both magnetic and resistive measurements. 
The crystals were plate-like, golden in color, and mea- 
sured ~ 500 x 300 x 20 /xm 3 in size. The crystal quality 
and orientation were characterized by Laue and indexing 
rotation photos, in which well-defined principal Bragg 
diffraction spots were identified and used to align the c*- 
and a*-axis to the scattering plane. Crystals chosen for 
the measurement showed rocking curve widths of ~ 0.05° 
for low order reflections. Beam damages to the crystals 
were not visible upon visual inspection using a micro- 
scope before and after the measurements. 

NIXS spectra were collected at room temperature on 
the Taiwan inelastic x-ray scattering beamline BL12XU 
at SPring-8 |l9| . The energy transfer was varied by scan- 
ning the incident energy relative to the near backscatter- 
ing energy (9.886 keV) of a 2-m radius Si(555) spherical 
analyzer [2(J . The total energy resolution was 65 and 250 
meV respectively using two configurations of the beam- 
line and the spectrometer (l9| . The momentum resolu- 
tion was 0.6 nm -1 in the horizontal scattering plane, and 
2.3 nm -1 in the vertical plane. 

Figure n shows a selection of the NIXS spectra taken 
over the energy region of the low-energy collective mode 
with momentum transfer q covering four entire BZ's 
along the c*-axis. All spectra are raw data normalized to 
the incident beam intensity in unit of counts per minute. 
The periodic nature of this mode in energy is apparent: 
The feature starts at ~ 3 eV at q = 3.0 nm -1 [Fig.QIa)], 
disperses down to ~ 4.5 eV at the first zone boundary 
at 8.9 nm -1 where the dispersion turns negative until it 
reaches near the next zone center at 17.8 nm , and then 
repeats itself in the third and the fourth BZ [Fig. IHb)]. 
The width of the loss feature varies also periodically, with 
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FIG. 1: NIXS spectra at various momentum transfers q || 
c*-axis showing the low-energy collective mode, where q — 
8.9 nm -1 corresponds to the first boundary of the extended 
Brillouin zones (BZ) . The total energy resolution was 65 meV 
for (a), and 250 meV for (b). 



about 0.5-eV (1-eV) width at half maximum near the 
zone centers (boundaries). Experimentally, the energy 
dispersion of this mode with momentum transfer can be 
described entirely by a simple cosine function (see Fig.0): 
lo = ujo — 27COs(gc), with u = 3.55 eV, 7 = 0.49 eV, and 
c = 0.352 nm the lattice constant of MgB2 along the o 
axis. An extrapolated value of 2.57 eV at q = can thus 
be obtained, which is in excellent agreement with recent 
optical studies [2l| and provides additional support for 
the experimental confirmation of this collective mode at 
low q. 

In order to understand the physics behind the data 
presented in Fig. ^ we have calculated the dynamical 
structure factor S(q,u>) based on time-dependent den- 
sity functional theory (TDDFT) 22], including fully the 
CLFE. Being a quantity directly proportional to the scat- 
tering cross section of NIXS |9j , the dynamical structure 
factor 5(q, u>) is related to the density-response function 
X by the fluctuation-dissipation theorem, 

S(q,w) = -2WTm XGq , Gq (q-G q , W ), (1) 

where G q is the unique vector of the reciprocal lattice 
which brings q into the first BZ. In TDDFT, \ obeys 
the formally exact integral equation \ = \s + Xs{ v + 
fxc)Xi where \S is the response function for Kohn-Sham 
electron-hole pairs, v the Coulomb interaction, and fxc 
the many-body kernel xs is evaluated for a ground 
state obtained in the local-density approximation (LDA) 
using the LAPW method [2^]. The many-body kernel 
fxn is evaluated for the adiabatic extension of the LDA 
|24|, defining the TDLDA response, which in view of our 
NIXS data, is an excellent approximation j2{|. The in- 
tegral equation is then solved for the response matrix 
XG,G'(q _ G q ,w) Hi|. The calculated £(q,u>), obtained 
using a damping parameter of 0.17 eV in the evaluation 
of xs, is compared with experiment in Figs. [21 and [21 

The calculated periodic energy dispersion of the collec- 
tive mode represented by the dashed line in Fig.[2]can be 
seen to agree almost perfectly with experiment through- 
out the four BZ's investigated. In Fig. 01 the calculated 
5(q, lo) is compared in detail with the NIXS spectra at 
a few selected q values of the second period. Here no 
attempt has been made to convert the experimental raw 
data onto absolute units due to insufficient data range 
to apply the f-sum rule The comparison is made 

therefore qualitatively with a constant scale factor be- 
tween theory and experiment for all the q values shown. 
Nevertheless, the spectral behavior of the observed col- 
lective mode is reproduced very well by the calculation 
despite an overall underestimate of the spectral weight 
at the higher energy side. 

Now, the density-response matrix can be ex- 
pressed in terms of the inverse of the dielectric 
matrix as xg,G'(<1 ~ G q ,u) = - G q + 

G) I [e(q — Gq,w)]e G , — Sg,G' (■ For the latter, we 
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FIG. 2: (Color) Theoretical S(n, to) calculated in the present 
work in false color log scale as a function of energy and mo- 
mentum transfer showing the cosine energy dispersion of the 
low-energy collective mode. Filled squares and triangles mark 
the energy positions obtained from the NIXS spectra shown 
in Fig. whereas filled circles are data from another set of 
spectra taken with a total energy resolution of 250 meV. 
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FIG. 3: Qualitative comparison between the calculated 
S(q,u) and the NIXS spectra at selected q values of the sec- 
ond period. 



have derived the following exact result |28j . 

1 



[e(q- G,^)]^ G = 



eG q ,G q (q- G q ,w) 
+F(q,q-G q ;^)[e(q-G q ,c)] 1 (2) 

for q's in higher BZ's 29]. Here F(q, q - G q ;w) = 
-W(q,q - G q ;w)/e Gq ,G q (q - G qi w), where VF(q,q - 
G q ;w) = (Det [e] - e Gqi G q M Gq ,G q ) /M , , M 0)0 and 
^G q ,G q being the minor of the (G = 0, G' = 0) and 



(G = G q ,G' = G q ) elements of e, respectively. Here, 
the arguments of all quantities shown symbolically are 
(q- G qi w). 

Equation (JSJ brings out the physics of the charge re- 
sponse for q's in higher BZ's quite vividly. The first term 
corresponds to the response appropriate for homogeneous 
media; it gives rise to charge fluctuations involving inco- 
herent electron-hole pairs. The second term introduces 
all the CLFE. For q — G q near the zone center, collec- 
tive charge fluctuations (plasmons) are built into the re- 
sponse function [e(q — G q , uj)] q ; whether such excita- 
tion is actually realized in 5(q, uS) (via Eq. (JTJ) depends 
on the impact of the material- dependent coupling func- 
tion F(q, q — G q , u). 

In Fig. 01 the contributions to the 5(q, oj) of MgB2 
from the first and second terms of Eq. (J2J are shown 
for q = 32.7 nm -1 in the fourth BZ along the c*-axis. 
Clearly, in the energy region relevant for the collective 
mode associated with the single-particle excitation gap 
of the B tt bands [l3L the first term amounts to 
a weak background, while the second term features a 
sharp peak, which is typical of q — G q near the zone 
center (|q — G q | = 2.9 nm" 1 in Fig. 0J). Now as shown 
in the insert of Fig. 01 ReF is positive and nearly constant 
for energies around the collective mode, while InxF goes 
through a local zero. Such an F ensures that the sec- 
ond term of Eq. (J2J) is controlled by Im [e(q — G q , w)]~J,. 
In other words, in MgB2 F feeds the small- (q — G q ) 
physics into the charge fluctuations for large q's. The 
plasmon can be viewed as acting - via F - as a source 
(additional to the external field of the x-ray photon) driv- 
ing a collective charge fluctuation generated in a NIXS 
event for large q's. It is important to notice that this 
process is intrinsically periodic. Indeed, the fact that for 
higher BZ's the loss feature retraces its energy dispersion 
in the first BZ reflects the impact on 5(q, oj), via the cou- 
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pling function F in Eq. J3J, of the "periodic" condition 
eG.G'C 4 , where k is in the first 



£G.G'(q - G q 

BZ. 

It is also interesting to note that there is a non-trivial 
difference between the response in the higher BZ's and 
that in the first one. In higher BZ's the collective mode 
arises from the "propagation" of the excitation gap of the 
B 7r bands built into eo,o via the inverse dielectric matrix 
of the second term in Eq. (0. This is the essence of the 
CLFE, in this context. By contrast, in the first BZ the 
excitation gap is directly built into the scalar response 
l/eo,o(q,w) 113- 
Additional insight into the charge response of MgB2 is 
achieved in light of a simplified version of Eq. (J2J, cor- 
responding to the 2x2 dielectric-matrix model invoked 
by Sturm, Schiilke, and Schmitz [3(J in their study of 
the S(q, uj) of Si. For Si, F is complex with a negative 
real part; thus the CLFE yield only a weak Fano-type 
resonance, which highlights the novelty of the charge 
response of MgB 2 . As we show elsewhere j2^|, a 2x2 
dielectric-matrix model represents a fair approximation 
to the full calculations we have performed for MgB2 using 
dielectric matrices of rank up to 23 and an 18x18x18 k- 
mesh. (For q— G q near each higher BZ boundary a 3x3 
model becomes necessary [28|.) In such a 2x2 model, 
the response function from the second term of Eq. 
becomes simply l/eo,o(q — G q ,a->), which makes obvious 
the identification of the plasmon discussed above with 
the collective mode discussed in Refs. O and Hi for q's 
near the first BZ T point. Most importantly, the fact that 
a dielectric matrix of rank two (or three) describes the re- 
sponse of MgB2 rather accurately is traced to the layered 
structure of MgB2 , which yields a first shell of non-zero 
G's containing only two vectors, and to the fact that the 
length scale of the polarization in real space is determined 
by orbitals which are rather extended, leading to a suffi- 
ciently fast decay of £g,g'; therefore, increasing the rank 
of the e-matrix beyond two or three does not affect the 
physics, in contrast to the case of, e.g., transition-metal 
oxides O. 

In conclusion, the electron-hole degrees of freedom in 
MgB2 lead to a novel charge response: long-lived collec- 
tive charge-density fluctuations can be excited in a NIXS 
event involving large momenta - corresponding to length 
scales which are comparable with the size of the B p or- 
bitals. This is qualitatively different from the case of 
simple metals and covalent semiconductors discussed in 
Ref. [33 , the difference stemming from the fact that the 
coupling function F is ~ 60 times larger in MgB2. The 
CLFE were shown to dominate the physics of the charge 
excitations for large momenta, leading to the striking 
periodicity exhibited by the observed low-energy excita- 
tion as the momentum transfer steps through successive 
BZ's. Overall, the impact of the CLFE on the charge 
response of MgB 2 rivals their importance in materials 
involving confined geometries such as carbon nanotubes 



[3l| and superlattices |32j . Ultimately, the nature of the 
measured charge excitations stems from the layered elec- 
tronic structure of MgB 2 and the delocalized nature of 
the orbitals involved in the screening. Analogous physics 
should be at play in other layered compounds of current 
interest. 

Stimulating discussions with B. Friedman, B. C. Lar- 
son, J. Z. Tischler, A. Shukla, M. Calandra, and K.- 
D. Tsuei are acknowledged. A.G.E. acknowledges sup- 
port from NSF ITR DMR-0219332. Research at ORNL 
is sponsored by the DOE, Office of Science, DMS un- 
der contract with UT-Battelle, LLC. W.K. and A.G.E. 
acknowledge collaboration via the DOE/CMSN. This 
work was carried out at SPring-8 under experiments 
No. C03A12XU-1500N and C03B12XU-1503N, and was 
partly supported by the National Science Council of Tai- 



[1 

[2 
[3 

[4 

[5. 

ir 

[8 
[9 
[10 

[11 

[12 
[13 

[14 

[15; 



[16] 
[17] 
[18] 
[19] 



[20] 
[21] 
[22] 



Electronic address: cai@nsrrc.org.tw 

Present address: HP-CAT, Advanced Photon Source, Ar- 

gonne National Laboratory, Argonne, IL 60439, USA 

J. Nagamatsu et al, Nature (London) 410, 63 (2001). 

S. L. Bud'ko et al, Phys. Rev. Lett. 86, 1877 (2001). 

J. M. An and W. E. Pickett, Phys. Rev. Lett. 86, 4366 

(2001). 

J. Kortus et al, Phys. Rev. Lett. 86, 4656 (2001). 

A. Y. Liu, I. I. Mazin, and J. Kortus, Phys. Rev. Lett. 

87, 087005 (2001). 

H. J. Choi et al, Phys. Rev. B 66, 020513(R) (2002), 
and Nature (London) 418, 758 (2002). 
S. Tsuda et al, Phys. Rev. Lett. 91, 127001 (2003). 
W. L. McMillan, Phys. Rev. 167, 331 (1968). 
W. Schiilke, J. Phys.: Condens. Matter 13, 7557 (2001). 
H. Raether, Excitation of Plasmons and Interband Tran- 
sitions by Electrons (Springer, Berlin, 1980). 

C. -Y. Moon, Y.-H. Kim and K. J. Chang, Phys. Rev. B 
70, 104522 (2004). 

A. Floris et al, Phys. Rev. Lett. 94, 037004 (2005). 
W. Ku, W. E. Pickett, R. T. Scalettar, and A. G. Eguiluz, 
Phys. Rev. Lett. 88, 057001 (2002). 
V. P. Zhukov, V. M. Silkin, E. V. Chulkov, and P. M. 
Echenique, Phys. Rev. B 64, 180507(R) (2001). 
A hint of this periodicity was contained in the NIXS 
data of Galambosi et al. [l6| who reported spectra for 
the first two BZ's. Similarly, the work of Goldoni et al. 
fl^l included a suggestion of this small-q collective mode. 
However, the rich physics behind the large-g NIXS data 
reported here was not anticipated by these works. 
S. Galambosi et al, Phys. Rev. B 71, 060504(R) (2005). 
A. Goldoni et al, Phys. Rev. B 66, 132503 (2002). 
Y. Takano et al, Appl. Phys. Lett. 78, 2914 (2001). 
Y. Q. Cai et al, in Synchrotron Radiation Instrumenta- 
tion: Eighth Internationl Conference on Synchrotron Ra- 
diation Instrumentation, AIP Conf. Proc. No. 705 (AIP, 
New York, 2004), p. 340. 

D. J. Wang et al, in Ref.0, p. 869. 

V. Guritanu et al, Phys. Rev. B 73, 104509 (2006). 
M. Petersilka, U. J. Gossmann, and E. K. U. Gross, Phys. 



5 



Rev. Lett. 76, 1212 (1996). 

[23] P. Blaha et ai, WIEN2k (Karlheinz Schwarz, Techn. Uni- 
versitat Wien, Austria), 2001. 

[24] E. K. U. Gross, J. F. Dobson, and M. Petersilka, in Den- 
sity Functional Theory II, edited by R. F. Nalewajski 
(Springer, Berlin, 1996), p. 81. 

[25] As it turns out, in MgB2 the many-body kernel plays 
a small role; in fact, even ignoring it altogether yields 
spectra which are in good agreement with our IXS data. 

[26] I. G. Gurtubay et ai, Phys. Rev. B 70, 201201(R) (2004). 

[27] J. Z. Tischler et al, Phys. Stat. Sol. (b) 237, 280 (2003). 



[28] O. D. Restrepo and A. G. Eguiluz, to be published. 
[29] For q's in the first BZ, Eq. © becomes: [e(q, lu)]~* = 

l/eo,o(q,^)+F(q,q-G 2 ;w) [e(q, w)]q* G2 , where G 2 = 
(27r/c)(0, 0, —1) |2q|. The first term here corresponds to 
the response function considered in Refs. [l3| and [l4T| . 

[30] K. Sturm, W. Schiilke, and J. R. Schmitz, Phys. Rev. 
Lett. 68, 228 (1992). 

[31] A. G. Marinopoulos, L. Reining, A. Rubio, and N. Vast, 
Phys. Rev. Lett. 91, 046402 (2003). 

[32] S. Botti et ai, Phys. Rev. B 70, 045301 (2004). 



